Constraints on CME Evolution from in situ Observations of Ionic 

Charge States 

Jacob R. Gruesbeck, Susan T. Lepri, Thomas H. Zurbuchen 
Department of Atmospheric, Oceanic, and Space Science, University of Michigan 

Ann Arbor, MI 48109, USA 

j agruesb@umich . edu 

and 

Spiro K. Antiochos 
NASA Goddard Space Flight Center 
Greenbelt, MD 20771, USA 

Received ; accepted 



ABSTRACT 


We present a novel procedure for deriving the physical properties of Coronal 
Mass Ejections (CMES) in the corona. Our methodology uses in situ measure- 
ments of ionic charge states of C, O, Si, and Fe in the heliosphere and interprets 
them in the context of a model for the early evolution of ICME plasma, between 
2-5 R s . We find that the data can be fit only by an evolution that consists of an 
initial heating of the plasma, followed by an expansion that ultimately results in 
cooling. The heating profile is consistent with a compression of coronal plasma 
due to flare reconnection jets and an expansion cooling due to the ejection, as ex- 
pected from the standard CME/flare model. The observed frozen-in ionic charge 
states reflect this time-history and, therefore, provide important constraints for 
the heating and expansion time-scales, as well as the maximum temperature the 
CME plasma is heated to during its eruption. Furthermore, our analysis places 
severe limits on the possible density of CME plasma in the corona. We discuss 
the implications of our results for CME models and for future analysis of ICME 
plasma composition. 

Subject headings: atomic processes - plasmas - solar wind - Sun: coronal mass 
ejections (CMEs) 



1. Introduction 


Coronal Mass Ejections (CMEs) are an important source of plasma and magnetic 
field in the heliosphere and the drivers of space weather events near Earth. They are 
characterized by the rapid and explosive release of large quantities of coronal material 
entrapped in its magnetic field. The rate of occurrence of CMEs is directly tied to the 
solar cycle (Gopalswamy et al. 2004) making them one important consequence of solar 
activity. Ejections from the Sun are more prevalent during solar maximum conditions and 
occur nearly 10 times less frequently (Gopalswamy et al. 2004) in the deepest of the solar 
minimum. 

CMEs are observed remotely in many wavelengths including white-light observations, 
by the Large Angle Spectroscopic Coronagraph (LASCO) (Brueckner et al. 1995), for 
example. Indeed, many CMEs are observed to be closely associated with enhanced plasma 
heating and energization, often relating to flare events during and immediately following 
the eruption (Harrison 1995; Vrsnak et al. 2005). Extreme ultraviolet (EUV) and soft 
X-ray emissions are most often used to analyze such heating (Krucker et al. 2010). These 
energization events can release energetic particles which are observed in remote observations 
of soft X-ray (Dauphin et al. 2006), hard X-ray (Temmer et al. 2010; Goff et al. 2005), and 
gamma-ray (Lin et al. 2003). Energetic particles can further interact with plasma near the 
Sun. in particular near the chromosphere, and lead to chromospheric evaporation events. 
We can directly observe the heliospheric manifestations of the erupted plasma, the so-called 
Interplanetary CMEs (ICMEs) which are characterized by specific plasma characteristics 
(Zurbuchen & Richardson 2006). 

ICMEs are usually associated with magnetically dominated plasma, which are 
approximately in a Taylor state with smooth and near force-free magnetic fields (Burlaga 
et al. 1981). ICMEs observed near 1 AU often have bulk speeds that exceed that of 



the surrounding solar wind and therefore are associated with interaction signatures, 
such as compressions or shocks (Zurbuchen & Richardson 2006, and references therein). 
Additionally, there are important expansion signatures, such as steadily decreasing velocity 
profiles, and depressed proton temperatures and densities (Gosling et al. 1973; Richardson & 
Cane 1995). Unusual charge states of heavy ions in the ICME plasma is another important 
signature (Bame et al. 1979; Henke et al. 1998: Lepri et al. 2001; Lepri & Zurbuchen 2004: 
Rodriguez et al. 2004; Zurbuchen et al. 2003). In fact, these compositional signatures were 
recently shown to be more prevalent than initially assumed, making them a powerful tool 
in ICME identification (Richardson & Cane 2004, 2010). 

EDITOR: PLACE FIGURE 1 HERE. 

Figure 1, shows the plasma, magnetic field and compositional signatures of an example 
ICME as observed with the Advanced Composition Explorer (ACE) (Stone et al. 1998). 
The steadily decreasing velocity profile, the low density and depressed temperatures are all 
signatures of an expanding magnetic cloud. The ICME is magnetically dominated and has 
the key characteristics of a flux-rope configuration. 

Figure 1G-K show measurements of the ionic distribution of C, O, Si and Fe. Unusually 
high charge states in C. O, Si. and Fe are the most important compositional signatures of 
this ICIUE. and many others; indicative of enhanced collisional ionization near the Sun. In 
fact, more than 80% of all ICMEs identified exhibit elevated charge states in Fe (Lepri et 
al. 2001) and generally, exhibit elevated states in the other species as well. The oxygen 
charge distribution, which usually peaks at 0 6+ in the solar wind, peaks at 0‘~, which is 
also observed in the Q t+ /0 6+ ratio, and has important contributions of 0 8+ . Similarly, 
carbon charge distributions in ICMEs are observed to be strongly dominated by C 6 “, the 
fully ionized state. The most unusual ionic charge state related to ICME plasma is the 



5 

Ne-like charge state of Fe (q — 16)elearly visible in Figure 1H (Lepri et al. 2001). Similarly, 
we will later discuss that silicon tends to be observed preferentially with a He-like charge 
state, ionized to Si V2+ , as shown in Figure II. 

The exact configuration of the observed ionization state is affected by collisional 
ionization, from electrons colliding with the ions in the plasma, and excitation autoionization 
(Cowan & Mann 1979; Hundhausen et al. 1968; Hundhausen 1972), where electrons from 
the outermost shells are lost due to the relaxation of an inner shell electron. The observed 
ionization states are an important indicator of the CME’s thermal environment close to 
the Sun. They are formed through a freeze-in process within 1-5 solar radii (Barne et al. 
1974; Buergi &: Geiss 1986; Hundhausen et al. 1968; Hundhausen 1972). Very close to the 
Sun, while time-scales for recombination and ionization are much shorter than the plasma’s 
expansion time-scale, the ionization states of heavies are approximately in ionization 
equilibrium, reflecting the local electron temperature. Due to the CME’s expansion into the 
heliosphere, the electron density decreases rapidly and thus the recombination time-scale 
increases as well. At a given heliocentric distance - depending on the specific characteristics 
of every ion, and the electron environment - each ionic charge states freezes in and no longer 
changes. At this point, the charge-state is said to be ’’frozen in.” Due to the characteristic 
differences of atomic physics for each ion, the heliocentric distance where freeze-in occurs 
varies widely, and simultaneously observed charge-state distributions become a sensitive 
measure for the thermal evolution of the CME plasma near the Sun during the entire radial 
range where freeze-in occurs. 

The unique nature of ICME charge states is further illustrated by Figure 2, showing 
a comparison of a time-period of solar wind composition observed from 0000 UT to 0200 
UT on July 30. 2003 with an ICME plasma observation from 0200 UT to 0400 UT on 
January 27, 2003. The characteristic differences are representative, and can be observed 



in Figures 1 and 3. The solar wind observation is comparable to the descriptions in von 
Steiger et al. (1997, 2000) for typical slow solar wind. Since CMEs originate from inside 
a helmet streamer, it is more appropriate to compare their composition with that of slow 
wind rather than the fast wind that originates from deep inside coronal holes. 

It is evident from Figure 2, that there is a noticeable increase in ionic charge states for 
each element in the ICME plasma compared to the typical wind. There is also a significant 
change in the qualitative nature of many distributions. Specifically, a significant fraction of 
C becomes fully ionized; 0 shows large deviations from 0 6+ , a He-like charge state that 
dominates solar wind measurements; Si charge states are substantially increased for the 
ICME period and peak at Si l2+ , the He-like charge state; Fe transitions into a bi- modal 
charge-distribution peaked around Fe 10+ and Fe 16+ , the Ne-like charge state discussed 
earlier. 

EDITOR: PLACE FIGURE 2 HERE. 

The purpose of this paper is two-fold. First, we survey ionic charge states in ICME 
plasma to determine the prevalence of bi-modal heavy ion distributions. Secondly, we 
develop a model to determine physical constraints for the ICME expansion especially 
focused on bi-modal Fe charge state distributions, while simultaneously, qualitatively 
re-creating the C, 0, and Si observations. Using this methodology, we will develop 
constraints for the temperature and density evolution of CMEs near the Sun. 


2. Observations 

This research is enabled, in part, by plasma observations from composition sensors 
such as the Solar Wind Ion Composition Spectrometer (SWIGS) on board the Advanced 
Composition Explorer (ACE) (Gloeckler et al. 1998), as well as its sister instrument on 



Ulysses (Gloeckler et al. 1992). These instruments independently determine speed, mass 
and charge of incident ions and therefore enable the measurement of full charge distributions 
in ICMEs. The time-resolutions of these measurements are limited by counting statistics of 
incident ions, and are typically limited to 2 hours for ACE-SWICS. The 0 7+ /0 n+ ratio is 
typically measured with 1 hour time-resolution. For each time-period, ions are accumulated 
and analyzed largely using the methodology described by von Steiger et a.l. (2000). The 
analysis of Fe charge states is rather straightforward because its mass of approximately 56 
axnu substantially exceeds the mass of other heavy ions of similar abundance in the solar 
wind. C and 0 are the most abundant heavy ions in the wind and are therefore easily 
discernible. However, the Si measurements currently are limited to only charge states 
between 8+ and 124- due to significant overlap with other neighboring ions in measurement 
space. Charge-state calculations are unsuccessful approximately 1% of the time due to low 
counting statistics or operational events on ACE independent of the solar wind plasma and 
compositional properties. 

The analysis of each ICME includes composition data, as well as plasma and field 
observations, as demonstrated in Figure 1. The proton moments obtained by the Solar 
Wind Electron, Proton and Alpha Monitor (SWEPAM) (McComas et al. 1998) data are 
averaged over a 64 second time period. The magnetic field observations are provided by the 
Magnetic Field Experiment (MAG) (Smith et al. 1998), averaged over 64 second intervals 
as well. ICME times and plasma boundaries are determined and given from Richardson & 
Cane (2010), which were determined primarily by magnetic and plasma signatures, while 
using composition and charge state data to confirm identification. 

Richardson k Cane (2004) have previously discussed the observations of generally 
enhanced ionization of C and O within ICMEs as compared to solar wind at the same 
speed. Indeed, such signatures were observed in 70% - 90% of all CMEs irrespective of 



their field configuration. ICMEs associated with a magnetic cloud are shown to have a 
higher concurrence of compositional anomalies than non-cloud ICMEs. 

Here, we focus on the bi-modal nature of Fe distribution, generally with peaks at Fe l0+ 
and Fe 16+ . Bi-modal Fe charge states are often observed in ICMEs, but they don’t always 
extend throughout the entire duration of the ICME. For example, bi-modal distributions 
are observed in Figure 1 during the first 4 hours, from time 0400 UT to time 0800 UT. The 
bi-modal nature of the ICME is not easily discernable from time 0800 UT till the beginning 
of February 19, but is seen again from then on until the end of the ICME passage on time 
1600 UT. Yet, bi-modal Fe charge states are equally likely to be observed during any time 
periods within ICMEs. Figure 3 shows ACE observations of a second ICME, observed on 
January 27 2003, presented in a format identical to that of Figure 1. The ICME boundaries 
were once again taken from Richardson k Cane (2010). 

EDITOR: PLACE FIGURE 3 HERE. 

The plasma composition during this ICME looks similar to the previously presented 
ICME data, in Figure 1: All four ionic distributions presented indicate significant 
enhancement during the event leading to higher average charge states, indicative of 
temperature enhancements. Also, the Fe charge state distribution is bi-modal for a majority 
of the ICME’s duration. 

It is well-known that ICMEs often include identifiable sub-parts (Mulligan & Russell 
2001). with each having specific magnetic field configurations or plasma characteristics with 
well-defined transitions between them. The ICME in Figure 3 includes such a qualitative 
transition at time 1600 UT on January 27. Prior to this, the magnetic field configuration is 
that of a magnetic flux-rope with the characteristic suppressed magnetic field turbulence, 
and reduced density characteristic of a magnetic cloud. The plasma temperature exhibits 



some interesting sub-structure during this time-period. At 1600 UT we observe a distinct 
transition visible in all plasma quantities - except for the plasma velocity. After time 1700 
UT. the magnetic field is again rather smooth, but following an orientation different than 
the one before. 

The compositional signatures appear to reflect this transition as well. SWIGS does not 
have a high enough time-resolution to analyze the transition with the same time-accuracy as 
SWEPAM and MAG, however there are discernible differences between the compositional 
signatures before and after time 1600 UT. The heavy ion composition during the first 
part of the ICME has substructures which have qualitative correspondence with variability 
in plasma temperature. It has also been argued that this plasma heating is strongly 
associated with flares (Lepri & Zurbuchen 2004; Rakowski et, al. 2007; Reinard 2005). 
Based on this observation, as well as analogous observations in Figure 1 and in many other 
ICMEs studied, we conclude that compositional structure reflects the plasma-structure and 
topology in ejecta. 

So. how prevalent are bi-modal Fe distributions in ICME plasmas? We address 
this question by surveying all ICMEs from 2000 till 2007 identified and characterized by 
Richardson & Cane (2010). For each ICME in this time-period, the percentage of bi-modal 
plasma is determined. This percentage is defined as the ratio of the time with bi-modal Fe 
observations divided by the total time during the ICME period during which composition 
data exist. The bi-modal characteristic is assumed to be present if there was a significant 
minimum in the ionic charge distribution, which includes more than one single charge state. 
This limits spurious minima caused by poor statistics in the data. 


EDITOR: PLACE FIGURE 4 HERE. 


Figure 4 shows the results of this survey, including approximately 200 ICMEs. Figure 
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4B shows a probability distribution of ICMEs according to the fraction that is observed in 
the bi-modal state normalized to 1. Figure 4A shows the cumulative probability distribution 
function (CDF) computed from 4B. It is immediately noticeable that more than 95% of 
ICMEs in this study include a measurable fraction of bi-modal Fe plasma. Using the CDF 
plot in the top panel, it can also be seen that 50% of ICMEs have at least 55% bi-modal Fe 
charge distributions. This bi-modal condition is thus found to be common in most ICMEs 
and therefore becomes a critical observational criterion for our understanding of thermal 
properties of CMEs near the Sun. 

It has been previously argued that bi-modal Fe observations are the result of plasma 
mixing between hot and cold populations (Gibson et al. 2006; Gopalswamy et al. 2001; 
Rakowski et al. 2007). However, our analysis together with the qualitative analysis of 
transitions in ICMEs, suggest that bi-modal charge distributions are an intrinsic property 
of a given plasma flow and not likely a mixing signature, which might be more diffusive in 
nature and more focused on specific boundaries and locales between cold and hot parts of a 
given CME. 


3. Model 


A freeze-in code has been developed, which solves for the final charge state distribution 

of a given atomic species under simple assumptions on a given CME’s heating and expansion 

z 

properties. To obtain the relative density of a given ion, y t = n,/ y) rij . and in particular 

z ‘=0 

the freeze-in condition, the following continuity equation, valid in the rest-frame of the 
ejected CME (Ko et al. 1997), is solved. 


dn, 

1 W 


llfil/i- I ( ' i- ! 


(T e ) - yi (Ci(T t ) + R^m) + y i+l Ri{T e )) 


( 1 ) 
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It is important to note that ra* is the number density of charge state i of the atomic 
species of interest, while y, is the relative number density. To solve this equation, the 
ion-specific recombination rates. R t , and ionization rates, C 2 , from Mazzotta et al. (1998) 
are used. These rates include electron ionization and auto-ionization as well as radiative 
and dielectronic recombination. These depend only on the electron temperature T e and are 
valid for all temperatures greater than 10 4 K. Photo-ionization is neglected here, as it plays 
a significant role for only singly and some double charged ions (Marsch et al. 1995). 

Equation 1 can be solved using a given electron density and temperature history (n e (t), 
T e (t)), which can be assumed, or constrained by observations. The system of ordinary 
differential equations described by Equation l is generally stiff as recombination and 
ionization rates can vary over many orders of magnitude. To tackle this, we use a 4 th 
order Runge-Kutta method that is specially suited for stiff equations (Press et al. 2002). 
The method also uses an adaptive step-size in order to use computational resources most 
efficiently. At each time, the method solves for the number density of an atomic species 
ionized to charge state i. We also check the accuracy of the integrator by testing the 
condition ]T)yj = 1. Our integrations are accurate to better than 10~ 6 . 

To model the heating and expansion of the plasma, we calculate using the equations 
above, the CME bulk parameters at each time t. We examine two types of evolution, one 
in which the CME plasma close to the corona is assumed to be in an elevated temperature 
and density state and then undergoes a pure adiabatic expansion. We show below that this 
evolution is not compatible with the composition data. The other evolution is one in which 
we include an initial rapid heating phase while the plasma is in the lower corona, and then 
impose the adiabatic expansion. Such a heating would be expected from the energy release 
due to the flare reconnection- that forms the ejected CME flux rope in the standard model 
(e.g.. Lynch et al. 2004). 
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Since the heating of the ejected plasma is not well-known from first-principles models, 
we assume a simple, as hoc form for this heating that has the minimum possible parameters. 
The plasma is taken to have some initial plasma temperature To at t = 0, which rapidly 
increases to a set maximum temperature (T max ), in a set heliospheric distance (r heat ). 
Generally, the heating distance is chosen to be some value between 0-0.5 R Sun from the 
surface. In addition to the heating profile, some assumption must be made as to the 
evolution of the density. The simplest is to take the density to evolve adiabatically, in 
other words, the heating is due to a pure compression of the plasma. Such a compression 
is, indeed, seen in simulations of CME/flare evolution. The reconnection jets coming out of 
the flare current sheet drive plasma compression both in the flare loop system remaining 
in the corona and in the escaping CME plasmoid. However, we also expect there to be 
some direct heating of the plasma due to magnetic dissipation throughout the CME/ICME 
evolution. This direct heating is highly model specific; therefore, to keep our results in 
this paper as general as possible, we consider only an adiabatic compression for the initial 
heating evolution, as well as an adiabatic expansion for the ejection. Of coure, given some 
quantitative model for CME temperature and density evolution, it would be straightforward 
to use it, instead of the adiabatic assumption, in our analysis. 

Equation 2 shows the forms of the temperature and density evolution that we use 
for the initial heating. Note that the model has only four important parameters: the 
initial temperature and density, the maximum temperature, and the radius at which this 
maximum temperature occurs. This is the minimum possible number of free parameters for 
any such model. 


T{t) 

n(t) 




T(t-) 

T max 


3/2 


where r(t) < r heat . 


( 2 ) 
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The particular functional form for the temperature evolution, a quarter sinusoid, is not 
critical. It was chosen merely to ensure that the temperature has a continuous derivative, 
which is necessary for greater accuracy of the numerical method. The evolution of the 
density is derived directly from the adiabatic expansion formula, assuming a monoatomic 
ideal gas. After the rapid heating phase, the plasma cools as its volume expands. Again, 
the adiabatic equations couple the density and temperature evolution. Equation 3 shows 
the equation which describes the evolution of these parameters. 



I where r(t) > r heat . 


( 3 ) 


Where f3 is an expansion factor of the plasma, typically set to a value between 2 and 
3. The cooling is dependent on the adiabatic index, 7, which is set 5/3, the value for an 
ideal monoatomic gas. Note that n rnax and T rnax occur at the end of the heating period and 
define the beginning of the expansion cooling phase. 

Motivated by observational characteristics of CMEs, our model allows for a non- 
constant velocity profile of the plasma. It. is particularly important to model rapid 
acceleration of CMEs near the Sun. CMEs rapidly accelerate to a maximum velocity, which 
remains nearly constant, much like the velocity profile of the solar wind (Gallagher et, al. 
2003). Using LASCO coronagraph images, a linear acceleration can be determined for 
some CMEs (St. Cyr et al. 2000). The SOHO/LASCO CME catalog makes these values 
available for CMEs when the calculation is possible. In our model, therefore, the CME 
plasma undergoes linear acceleration from the corona until the velocity reaches the observed 
in situ velocity, from ACE. After this, the velocity is set to the observed value and remains 
constant. 


Note that the total mass of CME plasma is assumed to be constant during the evolution: 



there is no net gain or loss of particles relative to the proton density. Furthermore, as 
the plasma is continually ionized, the free electrons remain in our plasma volume. This 
allows the plasma to experience as much recombination as made possible by the available 
electrons, creating quasi-local ionization and recombination. In addition, we have no specific 
requirements on the velocity distribution of the electrons. The model assumes all electrons 
are in the thermal population. Finally, we initialize the charge state of our plasma from the 
local thermodynamic equilibrium solution to the ionization equations discussed earlier. The 
model’s initial temperature is used to determine the ionization and recombination rates. 

To study the various physical characteristics of freeze-in within CMEs, we study a CME 
with an associated linear acceleration. This allows for a velocity profile to be determined. 
Specifically, we use a velocity profile matching the January 27, 2003 CME determined from 
the linear acceleration (St. Cyr et al. 2000) and ACE observations and whose velocity 
profile is shown in Figure 5. We then perform a parametric study to understand the 
temperature profiles that may lead to the observed characteristics. To demonstrate the 
impacts of the various model parameters on the freeze-in ionization rates, we perform a 
series of combinations of the following three characteristics which we found to be crucial 
for this calculation: rapid heating, rapid expansion and an elevated initial bulk density. 
We want to show that the absence of one of these key characteristics leads to qualitative 
disagreements to the observations summarized in Figure 2. The key assumptions of three 
illustrative calculations are shown in Figure 5. 

First is shown the case with rapid heating and rapid expansion, but with relatively low 
initial density, close to typical solar wind densities (Wilhelm et al. 1998). Second, results 
from the case having rapid expansion and an enhanced bulk density is shown. However, 
no heating is present and the plasma only cools. Finally, the results are shown when all 
three ingredients are included. This is the only case that reproduces, qualitatively the 
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observed freeze-in distributions. To determine the relative accuracy of these cases, we use 
the January 27. 2003 event shown in Figures 3, but it should be noted that our model is 
not sophisticated enough to match the data in quantitative detail. We seek only to match 
the general characteristics of the observed ionic charges states of heavy ions, such as the Fe 
bi-modal distribution. 


EDITOR: PLACE FIGURE 5 HERE. 


Figure 5 summarizes the basic assumptions of the input profiles the model uses; 
different color lines represent a different case of the model. The red line is the case with high 
density, rapid heating and rapid expansion. The blue line represents the case with rapid 
heating and expansion, but a lower initial density. The cyan line represents the case where 
there is an enhanced initial density, but no rapid heating. The top panel shows the density 
evolution for each model. The second panel shows the temperature evolution. Here, the 
red and blue lines are on top of each other, since they both experience rapid heating, while 
the blue line diverges initially. The initial temperature of the cooling only case was chosen 
such that it becomes similar to the temperature curves of the other cases downstream of 
the corona. This was done to insure that any difference in results was due only to the near 
coronal heating profile and not to any other portions of the temperature profile. It is also 
important to note that the shape of each case’s temperature curve qualitatively matches 
the shape of its respective density curve. This is a result of the temperature and density 
evolutions being tied together by the adiabatic relation. Finally, the last panel shows the 
velocity, which is identical for all three cases, as this curve was determined from LASCO 
and ACE observations of this particular CME by the previously described method above. 
We show below the final charge state results and evolution for these three model eases. 
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4. Results and Discussion 

As discussed in the previous section, we ran our freeze-in and CME expansion model 
for three different cases attempting to recreate a charge state observation from the January 
27, 2003 CME, from ACE, as shown in Figure 3. The first case simulates CME plasma 
which undergoes rapid heating near the corona, along with rapid expansion of its volume. 
However, the initial electron bulk density is low, at 6 x 10 7 cm' 3 , which is in qualitative 
agreement with nominal solar wind densities at the corona (Wilhelm et al. 1998). Figure 
6 shows the resulting charge state distributions of four atomic species ( C, O, Si and Fe) 
plotted at specific time-intervals during the expansion. Expansion times were translated 
into radial scale using the integrated velocity profile in Figure 5. 

EDITOR: PLACE FIGURE 6 HERE. 

The bottom four panels show, of figure 6, the evolution of the charge state distribution, 
for each species, from the coronal surface to a distance of 10 solar radii, after which the 
ionic charge states remain frozen-in for all considered cases. In fact, charge states are 
generally frozen-in within a distance of 5 solar radii (Buergi & Geiss 1986; Hundhausen 
1972), and our model agrees with this result showing the charge states to be frozen-in at 
a distance of approximately 3 solar radii. After this point, the charge states distribution 
remains constant. The top series of panels show the final resulting and frozen-in charge 
state distributions observed at 1 AU, the approximate location of in situ composition 
observations. It is immediately apparent that the resulting charge states are much too cool 
as compared to the data in Figures 2 and discussed in Section 2. Without a significant 
density increase over solar wind profiles, the resulting ionic charge distributions are not 
sufficiently ionized, leaving the charge states of all the species at values much lower than 
what is observed in ICMEs. 
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Next, we present the resulting charge state distribution from the case with initial high 
density and initial high temperature which only models the rapid expansion of an ICME. 
Thus, the CME plasma is assumed to be initially at a very high temperature. Figure 7 
shows the results in a similar manner as Figure 6. 


EDITOR: PLACE FIGURE 7 HERE. 

It is evident that the resulting plasma in this case of pure expansion has more enhanced 
charge states for all atomic species than the previous low density case. The problem, 
however, is that the charge states at 1 All are much hotter than is observed. The C, 0. and 
Si are all essentially fully ionized. The reason is that the plasma has simply not had enough 
time to recombine down to the observed charge states due to the rapid expansion and 
accompanying density decrease. We note that the Fe distribution does exhibit bi-modality, 
but the peak at Fe 16+ is much stronger than that at Fe 10+ , which is not in agreement with 
the observations. We conclude that rapid expansion combined with high initial density to 
be critical ingredients for producing bimodal Fe distributions, but that these two ingredients 
by themselves cannot be made consistent with all the compositional data. 

Finally, we present the case where all three ingredients (high density, rapid coronal 
heating, and rapid expansion) are present in the model. These results are shown in Figure 
8, again in a similar manner as the two previous figures. 


EDITOR: PLACE FIGURE 8 HERE. 


Figure 8 also shows a freeze-in behavior comparable to Figure 7. Again, charge states 
freeze in within the first few solar radii from the corona. Again, the model manages to 
match the CME like states of C and O as well as key characteristics of Si. However, under 



18 

this set of assumptions, a qualitatively similar bi-modal Fe distribution results, creating 
a consistent set of results under the assumption of one single density, temperature and 
velocity profile. Specifically, enhanced C and O charge states are seen, yet neither are fully 
ionized. Additionally, we see enhanced silicon charge states, and two distinct peaks in iron 
at Fe +1() and Fe +lb . The formation of bi-modal iron peaks is in part a result of high initial 
electron density and rapid expansion along with rapid heating, as can be seen by comparing 
figures 7 and 8 above. 

Figure 9 compares the results in Figure 6 and 8 with representative observations of the 
January 27 2003 CME from 0100 to 1400 hours on January 28, which again represents the 
behavior already discussed. 

EDITOR: PLACE FIGURE 9 HERE. 

The top row of charge state distribution in Figure 9 shows the result of the model 
which had rapid heating and expansion, but lacked a high initial density, the results of the 
case shown in Figure 6. The bottom row, from Figure 8, shows the result when the density 
is high and rapid heating and expansion are present. We left out the case shown in Figure 
7 because peaks are generally in the same place, and only the relative number density in 
the peaks differ. In addition, the Fe distribution of figure 7 is primarily uni-modal. The 
middle row shows the actual charge state distribution observations from ACE/SWICS of 
the January 27. 2003 ICME. which we are modeling here. All the charge states in the 
top-line are much lower than what is observed, while the relative shapes of the results in 
the bottom row are in qualitative agreement with observations, perhaps to a lesser degree 
for Si. Silicon from our model has a single peak, around charge states Si +9,+l °, while ACE 
observes a bi-modal silicon distribution with peaks at Si +9 and Si +}2 . Our model results 
in a considerable amount of Si +U ; which would be ionized further if the temperature curve 
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was modified slightly to have a shorter heating period. Finally, experimental issues limiting 
Si data were already discussed previously and may also effect this comparison. 

In order to test the robustness of our model’s code, we verified that a population of 
hydrogen would be fully ionized, using the third cases assumptions. It was found that 
hydrogen became fully ionized, as expected. Additionally, we repeated this test with for 
helium. Again, we found that helium becomes fully ionized; however, a small population 
recombines in the He + population, 10~ 4 in relative density. This is most likely due to 
the numerical solver, as this small recombination only occurs once the temperature of the 
plasma drops below the range of valid temperatures for the reaction rates. As previously 
stated, the rates used in the model are valid only for temperatures greater than 10 4 * * * A" 
(Mazzotta et al. 1998), however, during the full evolution of our plasma, the temperature 
does fall below this value. 

We also conducted tests to put a limit on how much the maximum density can 
vary with the particular temperature profile we have chosen. For the third case, which 
successfully recreated a bi-modal Fe distribution while matching C and O distributions, 
an initial electron density of 6.0 x 10 9 [cm -3 ] was used. By incrementing this value up 
and down we were able to determine a range of valid densities resulting in the Fe bi-modal 
distribution, which is 4.75 x 10 9 < n, < 8.5 x 10 9 cm' 3 . 

5. Conclusions 

The pervasive nature of the observed charge state characteristics at 1 AU (Lepri et 

al. 2001; Richardson <k Cane 2010) suggest a simple and ubiquitous physical process or 

explanation. From running a large number of test models, we conclude that charge state 

distributions are powerful tools for constraining the thermal properties of CME plasma near 



20 

the Sun. We find that bi-modal Fe charge distributions can evolve from a single plasma 
that rapidly expands from a high initial density. We also find that bi-modal Fe observations 
and elevated Si observations can only be made consistent with simultaneous C and O data, 
if the elevated initial density goes through a rapid heating before it cools in the expansion. 
No further assumptions are needed to fit qualitatively these data. 

Of course, we cannot exclude the possibility that other candidate processes, such as 
non-thermal electron characteristics near the Sun, or more complex evolutions involving 
multiple sources of plasma can achieve similar agreement with observations. We would 
argue, however, that the basic plasma evolution derived in this work, a rapid compression 
followed by an expansion, is the simplest and most likely explanation for the observations. 
In order to observe Fe lfi at 1 AU, the plasma must reach fairly high temperatures, several 
MK, before the radius at which the freeze-in condition sets in. But if the plasma is in 
thermal equilibrium at this point, then the resulting charge states for the lighter elements 
would be too high at 1AU, as in figure 7. The only way to avoid this discrepancy is to heat 
the plasma rapidly so that the lighter elements never reach their fully ionized state, while 
the iron achieves high ionization. Furthermore, since the velocity and expansion properties 
of CMEs are fairly accurately known from coronagraph observations, the observed charge 
states distributions impose tight constraints on the plasma density during this heating and 
expansion evolution. Too low a density implies that the various elements never reach the 
observed ionization states; whereas too high washes out the bi-modal peak in iron. These 
straightforward and compelling arguments demonstrate the power of our compositional 
analysis techniques for deriving the detailed properties of the solar origins of space weather 
from in situ measurements in the heliosphere, 

Note that the evolution determined by our composition analysis: a rapid compressive 
heating to high density followed by expansion cooling is fully consistent with the expected 
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effects of flare energy release on the thermodynamics of CME plasma (Harrison 1995 : 
Vrsnak et al. 2005). These effects have been discussed by many authors (e.g.. Canfield & 
Reardon 1998; Li et al. 2008) and typically relate to reconnection processes going on near 
the Sun (Lynch et al. 2004). The upward directed reconnection jets from the flare current 
sheet are expected to compress and heat the CME plasma, as derived above. By combining 
our data and modeling techniques with sophisticated MHD simulations (Antioehos et al. 
1999; Gombosi et al. 2000; Linker et al. 1990), we can obtain a powerful tool for the analysis 
of CMEs and their evolution near the Sun. The computations presented in Figure 8 provide 
predictions of the presence of certain ions at specific times during the CME’s evolution and 
should be testable, especially by forward modeling in which the calculation of the charge 
state evolution is included in the MHD simulation. A key point is that the various models 
for CME onset predict different locations for the initial heating of the coronal plasma that 
is ejected into the heliosphere as an ICME (e.g., Forbes et al. 2006). Of course, many of 
these differences are expected to be washed out during the propagation to 1 AU; however, 
with the advent of composition measurements near the Sun from missions such as Solar 
Probe Plus, our analysis may be able to determine not only the coronal evolution of CME 
plasma, but the fundamental mechanism responsible for the eruption itself. 

This work used the SOHO/LASCO CME catalog developed by NASA, NRL and the 
Catholic University of America. SOHO is a project of international cooperation between 
ESA and NASA. This work was performed, in part, with support from NASA grants 
NNX08AI1 1G, NNX07AB99G, and NNX08AM64G, as well as the GSRP program. THZ 
acknowledges the hospitality of the International Space Science Institute, where much of 
this work was being performed. SKA was supported, in part by the NASA HTP and 
TR&T Programs. The authors would like to thank Eric Christian for helpful discussions 
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Fig. 1. Summary of plasma parameters and charge state distribution for an ICME observed 
on February 18. 2003 by the ACE spacecraft. From top to bottom we show proton velocity 
(Vp). proton density (N p ), proton temperature (T p ), magnetic field magnitude (>£>]). RTN 
longitude (A), RTN latitude (5), () 7 ~ /()" " ratio, Fe charge state distribution }. Si charge 
state distribution (Qsi), O charge state distribution (Qo), and Ca charge state distribution 
(■ Qca )■ The ICME plasma field begins at the solid red line, at 0400 UT on February 18, 
2003, and ends at the. dashed red line, at 1600 UT on February 19, 2003. It can be seen that 
the all four ion experience enhancement in the ICME plasma. 
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Fig. 2. A comparison of the normal solar wi 
ICME plasma shows enhanced charge states f 


CME Plasma 
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(left) and of ICME plasma (right). The 
all four atomic species compared to the 


typical solar wind. Note the bi-modal structure of Fe with peaks at Fe w+ and Fe 16+ . For 


more details, refer to text. 
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Fig. 3. Summary of plasma parameters and charge state distribution for an ICME observed 
on January 27. 2003 by the ACE spacecraft. The ICME plasma field begins at the solid red 
line, at 0100 UT on January 27, 2003, and ends at the dashed red line, at 1400 UT on 
January 28, 2003. All ionic charge distributions show a significant enhancement during the 
ICME plasma. It can also be seen that a majority of the Fe distribution displays a bi-modal 
shape. Xot.e the qualitative compositional change at time 1600 UT of .January 27, when the 
magnetic field abruptly changes, denoted by the gold line. A bi-modal characteristic can 
also be observed in the middle of day 28. 
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Fig. 4. Statistical survey of bi-modal Fe distributions in ICME plasma. The bottom panel 
is a histogram showing the results of the survey. Percentage of total CMEs investigated is 
the vertical axis, while percentage of bi-modal plasma in the CME is on the horizontal axis. 
The top panel shows a cumulative distribution function (CDF) of the resulting dataset. This 
vertical axis is the probability that a given CME will have at least percentage of bi-modal 
Fe in the plasma field as the corresponding percentage on the horizontal axis. 
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Charge State 

Fig. 6. The results from the model which includes rapid heating and expansion, but with 
the omission of an enhanced initial density. The bottom set of panels show the evolution 
of the charge states from the coronal surface to a distance of 10 r gun away plotted for 
arbitrary instances during the expansion. For all four atomic species, it can be seen that 
freeze-in occurred around 3 r Sun . The top panels show the frozen-in charge state distribution 


observed at 1 AU. 
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Fig. 7.— The results from the model which includes rapid expansion and also has a suffi- 
ciently large initial bulk density, but omits rapid heating. Format as in Figure 6. 
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Fig. 8. The results from the model for the case where the plasma is rapidly heated, 
experiences rapid expansion, and has a sufficiently high initial density. Format as in Figure 
6 . 
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Fig, 9. Charge state distributions for carbon, oxygen, silicon, and iron. The top row is 
the final distribution result from the model in Figure 6. The middle row is the observed 
ACE/SWICS charge state distribution for a two hour span of the January 27. 2003. The 
final row is the final distribution from the model in Figure 8. 





